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Abstract 

We investigate a density-functional theory (DFT) approach for an unpolarized trapped dilute 
Fermi gas in the unitary limit . A reformulation of the recent work of T. Papenbrock [Phys. Rev. A, 
72, 041602 (R) (2005)] in the language of fractional exclusion statistics allows us to obtain an 
estimate of the universal factor, £3^, in three dimensions (3D), in addition to providing a systematic 
treatment of finite- N corrections. We show that in 3D, finite- N corrections lead to unphysical 
values for £3^, thereby suggesting that a simple DFT applied to a small number of particles may 
not be suitable in 3D. We then perform an analogous calculation for the two-dimensional (2D) 
system in the infinite-scattering length regime, and obtain a value of = 1- Owing to the unique 
properties of the Thomas- Fermi energy density- functional in 2D our result, in contrast to 3D, is 
exact and therefore requires no finite- A corrections. 
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I. INTRODUCTION 



Harmonically trapped ultracold atomic Fermi gases have received considerable attention 



within the last few years owing to the experimental realization of Fermi degeneracy 
DeMarco and Jin I], the pioneering work of the Thomas KTOup in the unitary regime \ 
and other experimental efforts including, for example, Refs.y,0,y,Q. Using the phenomenon 
of a Feshbach resonance, the two-body inter-particle interactions between the atoms can be 
continuously tuned according to the magnitude of the magnetic field across the resonance 
region. Of particular interest is the so-called unitary regime which occurs at the midpoint 
of this crossover, and is characterized by the divergence of the scattering length due to the 
existence of a zero-energy bound state for the two-body system j^J . In this regime, the only 
relevant length scale is set by the Fermi momentum, kj , and so the corresponding energy 
scale is the Fermi kinetic energy. Therefore, in three-dimensions (3D), the energy density 
should be proportional to that of a free Fermi gas 
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where is a dimensionless proportionality constant, and kf = (37r 2 p) 1 / 3 . In the uni t an 
limit, the system details do not contribute to the bulk thermodynamic properties 
and so is usually called the universal factor. Owing to the lack of a small expansion 
parameter, i.e., (akf), the exact determination of the universal factor C, 3 d is a challenging 
problem because the usual Green's function techniques for the many body system are com- 
pletely unreliable. Moreover, as the unitary Fermi gas may be relevant for the physics of 
neutron stars and high-T c superconductivity, the extraction of the universal factor is also an 
important task. 



To date, theoretical interest in the universa 



3D unitary Fermi gas 
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factor has been primarlily limited to the 
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13, 



181 ] . These theoretical approaches 



vary in their level of sophistication, ranging from simple applications of density-functional 
theory (DFT), to relativstic QED-like theories with hidden local symmetry, and quantum 
Monte Carlo (QMC) simulations. Unfortunately, the wide variety of theoretical approaches 
employed for the determination of £30 have also led to a broad range of possible values for 
£ 3 £>; namely, theory predicts £ 3 £> ~ 0.3 — 0.6, although Engelbrecht et al. have established an 



upper bound of £ 3 £> < 0.59 from a BCS treatment 



12j . The disagreement between theoretical 
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approaches not withstanding, it is generally accepted that the recent QMC simulations 
of Carlson et al. yield the most reliable estimate for the universal factor, which is 
£ 3D = 0.44 ±0.01. 

Initially it was difficult experimentally to obtain a value for £ 3 £>, because the precise 
location of the required Feshbach resonance was not known and a range of estimates 



£ 3 d ~ 0.3 - 0.8 



19| were thus obtained. With the position of the resonance now accu- 



rately determined, more recent experimental measurements of fen are in better agreement 



2l|. These results strongly 



and have yielded £ 3D = 0.51 ± 0.04 |20| and £ 3D = 0.46 ± 0.05 
suggest that the QMC estimate of ~ 0.44 is indeed very reliable. 

In this paper we determine the universal factor, £ for a trapped two-dimensional (2D) 
and 3D Fermi gas. In 3D, our formulation of the problem in the language of fractional 
exclusion statistics (FES) allows us not only to extract a value for £3^, but also provides a 
systematic treatment of finite- N corrections. Our determination of £ 2 d is likewise obtained 
within an exceedingly simple mathematical framework, which nevertheless leads to an exact 
result for the universal factor in two- dimensions. We are aware of no experimental results for 
the (quasi) 2D Fermi gas in, or near, the appropriate scattering regime and only two recent 



theoretical studies with the explicit aim of extracting both agreeing on £ 2 d = 1 |22l. 12c 
It has also been pointed out in passing by Tonini et al. |24j that for large scattering length 
in 2D, the effective interactions are weakly attractive, thereby yielding the non- interacting 
Fermi gas in the infinite scattering length limit this observation plays a critical role in 
our subsequent determination of 



II. TRAPPED FERMI GAS IN THE UNITARY REGIME 

In this section, we apply a simple DFT in order to determine the universal factor for an 
unpolarized trapped Fermi gas. We begin by reformulating the work of Papenbrock in 3D, 
with particular attention paid to finite-iV corrections, and what impact they have on the 
estimate of £3^3 . We also present an analogous analysis in 2D, and show that it yields an 
exact result for £ 2 d without requiring any finite- N corrections. 
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A. Three Dimensions 



As stated in the introduction, Carlson et al. have performed simulations for a uniform 
3D Fermi gas in the unitary regime. Dimensional arguments suggest that in this regime, the 
energy must be proportional to that of a free Fermi gas, 

E(N) = £3D£free(iV) , (2) 

where Ef ree (N) is the Thomas- Fermi (TF) energy of N noninteracting spin- 1/2 fermions. 
Surprisingly, the QMC simulations with iV = 10 — 40 fermions agree very well with Eq. (J2J) 
in spite of the fact that the exact quantum mechanical energies are known to differ from 
the corresponding TF energies due to shell and finite size effects. Motivated by this result, 



and again on dimensional grounds, Papenbrock 
from TF theory, 



3 



has suggested that a density functional 



4 D F H °(p)=^^cp 5 / 3 + ^Vp, (3) 

should be a good approximation to the exact density functional. In general, gradient cor- 
rections should be included, but Eq. (3) represents the simplest possible form for the energy 
density-functional. In the above, u is the trapping frequency of an isotropic harmonic oscil- 
lator (HO) potential, and c = ^(37r 2 ) 2 / 3 . It is straightforward to show that a minimization 
of this functional, subject to a fixed number of particles, leads to 

i / o \ 3/2 / 2 \ 3/2 

where I = (h/mu) 1 ^ 2 is the oscillator length. Integrating this energy functional up to the 
TF radius yields the total TF energy 

« HO = Imnis^. (5) 

The energy can be calculated exactly for the case N = 2, where there is a binding energy 
of hjj in the unitary limit such that total energy for two fermions in a harmonic trap 



E ex = 2fkv 25]. Comparison with the TF expression (which is expected to be reasonable 
for all even numbers of fermions) with N = 2 gives a value for the universal constant of 
£ 3D 0.54. Although it is not explicitly mentioned in Ref. H3, if one uses the exact two- 



r 



particle density 



in Eq. (3) rather than Eq. (4), and then performs the spatial integration, a much better 
value of ~ 0.48 is obtained. Thus, a very simple DFT already yields a value for £30 that 
is within 8 % of the generally accepted QMC result. It is then natural to ask what effect 
finite- N corrections will have on this simple result; afterall, TF DFT is generally only valid 
for N 1, and/or slowly spatially varying potentials, which is certainly not the case in the 
present analysis. In order to investigate the finite- N corrections, it turns out to be more 
convenient to reformulate the DFT of Ref. ^ in a slightly different manner. 



Bhaduri, et al. 



261 ] have argued that the uniform 3D interacting Fermi gas at unitarity 



can be described in terms of a non-interacting gas obeying fractional exclusion statistics 
(FES). This generalization of the standard Bose-Einstein and Fermi-Dirac particle statistics 



was developed by Haldane 



in the context of the quasiparticle excitations of the fractional 

. In these Haldane- Wu statistics the 



quantum Hall effect, and later generalized by Wu 
concept of the Pauli exclusion principle is extended such that the distribution function for 
a single particle state with energy e is given by 

1 
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where the temperature dependent function w is defined through 



w(l + w 
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In the limit a = this then reduces to the usual Bose distribution and for a 
Fermi-Dirac distribution. At zero temperature the distribution reduces to 
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The parameter a that defines the statistics gives a modified exclusion principle which states 
that the maximum occupancy of a single state is 1/a - hence fractional exclusion statistics. 

Let us now consider the density of states, D(e), for the 3D HO (including the spin 
degeneracy factor of 2) 

e 2 1 

D(e) = 7— - — — . (10) 

Following Bhaduri et al. 



26| |. we now utilize the Haldane- Wu statistics to write the total 



number of particles, N, in the system as 



N = - I D{e) de 
a 



(11) 



where, as above, 1/a replaces the usual zero temperature occupancy factor of unity for 
fermions. Rescaling all energies with respect to the oscillator energy hv gives the following 
expression for the chemical potential upon integration 

£-£-(aJV) = 0. (12) 
This equation has one physically acceptable root, viz., 

1 
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2 (12aJV + Vl44a 2 iV 2 - l) 1 / 

= (3ajV)1/3 + \ WnW ~ 152 WnW + " ' (13) 

to 0(iV~ 5 / 3 ) corrections. Similarly, the total energy of the system is given by 
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E = — D(e)e de 
a Jo 

-!(t-S)- (14) 

Rewriting Eq. (pHj) explicitly in terms of N results in the expression (to 0(N~ 2 ^ 3 ) corrections) 

4 + 8 + 32 384(3iV) 2 /3 + 

We can now make a formal connection with Ref. [US by identifying = a 2//3 . Retaining 
only the first term in Eq. (15), we then obtain 

E = i(3iV) 4/3 (6D) 1/2 , (16) 

which is identical to the result for the total TF energy given in Q]. For the N = 2 case, we 
have E e ^ = 2, and we obtain as a first approximation to a 

2 = a V3^_ w 2.725a 1 / 3 , (17) 

which immediately gives a = 32/81, or ~ 0.54, in complete agreement with the DFT 
approach of Papenbrock. If we keep the first leading order correction in Eq. (15), we obtain 
(again, for the N = 2 case) a second approximation to a 

fi 4/3 c2/3 

2 = a^ b — + a-^ b — , (18) 

which has no real roots. However, if we retain the first and third (i.e., a _1 /32) terms in 
Eq. (15), we obtain a much improved value of £3/3 ~ 0.49. Nevertheless, beyond the first and 



this latter approximation, we find that there is no real value of a which yields E = 2. Thus, 
a systematic treatment of the finite-iV corrections (sometimes referred to as the extended 
Thomas- Fermi theory (ETF)), with the goal of yielding an improved value for £ 3 £>, fails 
in 3D (for the N = 2 case at least). Of course this could also be an indication that the 
identification of the strongly interacting Fermi gas with a gas of non-interacting particles 
obeying FES is not valid (in 3D) for small particle numbers once finite- N corrections are 
considered. Bhaduri et al. |26j have however used this approach, even for small numbers 
of particles, to obtain thermodynamic properties of the unitary gas at finite-temperature 
which compare well with the Monte Carlo results 3|. It is also worth noting that for 
a = 1 (i.e., Fermi-Dirac statistics), the exact energy and chemical potential are E = 3 
and fj, — 2, respectively. In this case, the systematic corrections discussed above do yield 
improved values. For example, the first approximation gives E = 2.73 and /J = 1.82, while 
the second approximation already gives E = 3.14 and /x = 1.96. This said, the fact that the 
first approximation, and its augmented form (i.e., including the iV-independent term) yield 
values for the universal factor so close to the QMC result is perhaps fortuitous. Interestingly, 
the requirement that E be real leads to a lower bound of > 0.12, and the energy takes 
on a minimum value E ~ 2.45 for £ 3 £> w 0.3. 

B. Two Dimensions 

Turning now to the two-dimensional isotropic harmonic trap, the TF density functional 
reads (for the noninteracting Fermi gas) 

4 D F HO (p) = £p 2 + ^Vp. (19) 



At zero temperature, when filling M + 1 oscillator shells, Brack and van Zyl [30j have 
provided an exact, closed form solution for the zero-temperature particle density in 2D, 

P( r ) = — jr ]T(M - n + l)(-l) n L n ( —r 2 J e~ r , (20) 

71 n=0 ^ ' 

where L n (x) is a Laguerre polynomial. Inserting Eq. (|20j) into Eq. (|19|) . and integrating over 
all space yields the total TF energy 

E 2Dho = ^[2M 3 + 9M 2 + 13M + 6] . (21) 



The number of particles in M + 1 filled shells is given by N(M) = M 2 + 3M + 2. What 
is truly remarkable about Eq. (|21|) is that it is also the exact quantum mechanical energy 
of the system. In other words, the TF functional (which is exact only in the homogeneous 
limit), without gradient corrections yields the exact total energy of the inhomogeneous non- 
interacting system. 

Before applying the above result to the trapped 2D Fermi gas at unitarity, we must first 
make clear what is meant by "unitarity" in two-dimensions. If one considers the 2D case (see 
our comments for 3D in 8j) the s-wave scattering cross-section takes the form (h = 2m = 1) 
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with a 2 D now the 2D scattering length. The maximal value for the cross-section is then 
obtained when the logarithmic term vanishes. Following the equivalent argument to that in 



3D, the unitarity limit is therefore obtained when 



a 2 Dk 



1. In genera 



22 and 



jthis is not what is 



23 refer to the limit 



taken as the unitary limit in the literature. For example, Refs. 
a 2 o —> oo. Clearly, if we hold that unitarity occurs when o 2 d reaches is maximal value, 
the dimensional arguments leading to fact that the energy of the interacting gas must be 
proportional to the energy of the noninteracting gas are not applicable. In particular, the 
hierarchy: a ^> kj 1 ^> r , implicity used in obtaining Eq. (2) for 3D, is not sufficient in 2D 
as the range of the potential, Tq, cannot generally be taken to be zero. Indeed, unitarity in 
2D implies that a 2 o = kj 1 . Thus, if we wish to preserve an analogous notion of unitarity 
in 2D, we require that the 2D scattering length become very large such that a 2 r>kf 3> 1. 
Obviously, this does not correspond to a maximum in a 2 D, but the large (i.e., essentially 
infinite) scattering length limit is what is often regarded as the unitary regime in two- 



dimensions 



22, 



23). 



With the above qualification regarding unitarity in 2D in mind, we now consider the 
binding energy of a dimer in free space, E 241 ]. 
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where 7 = 0.57721 is Euler's constant, along with the coupling constant for s-wave 

scattering between two particles, 
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where G = 0.91596.... is Catalan's constant. In the limit ci2D — > oo, for a fixed density 
and range ro, go 0~ and Eq — > + . Thus, in the limit a2D^/ 3> 1, the trapped 2D 
gas corresponds to a weakly attractive Fermi gas with a zero energy bound state. In other 
words, the threshold for the appearance of the first zero-energy bound state corresponds 
to zero coupling. In this sense, the infinite scattering length limit in 2D corresponds to a 
noninteracting Fermi gas, which is in stark contrast to 3D, where the {a^okj) S> 1 regime 
requires a non-trivial treatment of interactions. 

We are now in a position to apply the simple DFT approach to the trapped 2D Fermi 
gas in this limit. Arguing as in Ref. [13J, the appropriate 2D TF density functional is given 
by 

4 D F H °(p)=6^p 2 + ^Vp. (25) 

Again, gradient corrections generally have to be introduced, but as we will show, they are 
not required for the 2D case. Minimizing with respect to p, subject to the constraint that 
the total number of particles is held fixed, yields a chemical potential 

/4 D F H ° = il'nN^h, (26) 

with TF density 

= ("-° " i B ™ v ) ' (27) 

Inserting Eq. (J27|) into the energy functional and integrating to the TF radius yields a total 
TF energy 

2S™o = \N^hw. (28) 

As we have already discussed, in the infinite scattering length regime the exact energy for 
two fermions in a harmonic trap is zero energy bound state. Hence the exact energy for a 
pair of fermions in a 2D harmonic trap is also E ex = 2hw. Comparison of the TF energy, 
viz., (|28|) with Eq. (21), which is exact in the infinite scattering length limit, gives a value 
for ^2D such that (31 1 

= 1 + w (29) 

For example, with N = 2, we obtain £ 2 £> = 9/8 = 1.125. Thus, while an analysis of finite- N 
corrections in 2D is not applicable here, because D(e) = e/(hu) 2 , Eq. (29) illustrates that 
employing the 2D TF density in Eq. (25) leads to (what turns out to be) an overestimate of 
£ 2 d, with &d — > 1 + as N — ► oo. 



We have, however already established that Eq. ([19)1 yields the exact quantum mechanical 
energy provided the exact density is used as input. Since the 2D gas in the infinite scattering 
length limit corresponds to the noninteracting system, if we use the exact iV-particle density, 
instead of the TF density, viz., Eq. (J27|) . and then perform the spatial integration of (|19j) . 
we immediately obtain £ 2 d — 1 identically. This is true for any N = M 2 + 3M + 2, and 
thus establishes that the simple 2D TF functional does indeed lead to a universal = 1, 
independent of the particle number, when the exact density is used. Recall that in the 3D 
case, the use of the exact two-particle density, Eq. (6), does lead to an improved value of 
£ 3D pa 0.48 as compared to that obtained using the TF density. However, as the 3D TF 
functional is not exact in the unitary regime, neither is the extracted universal factor. 



III. CONNECTION WITH FES IN THE 2D FERMI GAS. 

In a relatively recent paper, Srivastava et al. have considered a gas of N fermions interact- 
ing via a repulsive two-body zero-range (psuedo) potential within a mean-field approach [32] . 
Specifically, they showed that within the TF formalism, the zero temperature density func- 
tional (with a spin- degeneracy factor of 2 included) is given by 

411° = «|V + \mu^p , (30) 

where a = (1 + f^jr) is a dimensionless statistical parameter, and Mq is the zeroth moment of 
the two-body potential. The above equation implies that a gas of fermions in 2D, interacting 
by a repulsive two-body zero-range potential, appear as particles obeying ideal FES (at least 
within the TF framework). Notice that the form of Eq. (30) is identical to Eq. (25), but 
here, a is a statistical parameter unreliant upon being in the unitary or infinite scattering 
length regimes. In fact, the scaling of the kinetic energy functional here relies entirely on the 
details of the two-body potential, whereas in the unitary limit, the details of the potential are 
unimportant. The similarities between Eqs. (25) and (30), and the analgous expressions for 
p, fi and E, naturally lead to the question of what connection the Fermi gas in the infinite 
scattering length limit has with FES. Unfortunatley, the answer is quite anticlimactic in 
2D because it is only in the limit a — 1, that there exists any correspondence. This, 
unsurprisingly, leads immediately to = 1 since a = + yj$- 
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IV. CONCLUSIONS 



In this paper we have reformulated the simple DFT of Papenbrock in the language of 
fractional exclusion statistics to obtain an estimate of the universal factor £ 3D which is in 
good agreement with that obtained via QMC simulations, and experiments. We have shown, 
however, that trying to go beyond this estimate by using a systematic treatment of finite N 
corrections fails in 3D. 

In contrast, when we apply the DFT technique to the 2D trapped Fermi gas in the infinite 
scattering length limit (and use the exact density in the 2D TF energy density functional), we 
have shown that the exact quantum mechanical energy is obtained for any number of closed 
shells. The TF result thus obtained is therefore not limited to N 3> 1 as in the 3D case and 
we are able to establish that the universal factor for the harmonically trapped 2D Fermi gas is 
exactly £20 — 1- We therefore conclude that a simple DFT approach, as originally proposed 
by Papenbrock, is only strictly valid in 2D, where the TF energy functional, Eq. (19), has 
the remarkable property of being exact, without gradient corrections. 
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